Forecasting systems

ABSTRACT

A computer-implemented method of determining an approximated value of a parameter in a first domain is described. The parameter is dependent on one or more variables which vary in a second domain, and the parameter is determined by a function which relates sets of values of the one or more variables in the second domain to corresponding values in the first domain. The method is implemented on a computer system including a processor, and the method comprises: determining a plurality of anchor points in the second domain, wherein each anchor point comprises a set of values of the one or more variables in the second domain; evaluating, at each anchor point, the function to generate corresponding values of the parameter in the first domain; generating an approximation function to the function by fitting a series of orthogonal functions or an approximation to a series of orthogonal functions to the corresponding values of the parameter in the first domain; and using the approximation function to generate the approximated value of the parameter in the first domain.

TECHNICAL FIELD

The present invention concerns improved forecasting systems, in particular, but not exclusively to a system for accurately carrying out predictive functions using approximation techniques.

BACKGROUND

Predicting future values of parameters is desirable in many fields and endeavours, including molecular biology, oil and gas, structural engineering, financial modelling and weather forecasting. Typically, in these fields, the value of the parameters is dependent on a complex relationship between several underlying variable factors and sub-variables. Forecasting (i.e. determining future potential events with a probability associated with each of them) is useful as decisions and expectations can be made based on predictions. Also, such predictive analysis can result in useful feedback being provided to a control system to configure a system or device to accommodate the future state or events.

Predictive functions using approximation techniques are used in situations where a given parameter can only be determined by consideration of a group of functions where each function is determined based on multiple scenarios. The scale of computational effort required for this future behaviour simulation is so large, for example there may be millions of variable values requiring millions of approximation calculations, that even current state of the art computing platforms struggle to meet the desired level of accuracy and performance. For example, current forecasting systems may reduce the number of approximation calculations performed, thus reducing accuracy, or increase the amount of processing power used, thus increasing the electrical power and resource requirements of the system. The application of such predictive analysis techniques includes fields such as structural engineering, molecular biology weather forecasting, oil and gas and financial modelling.

For example, in biomedical and biological research, it is of the utmost importance to determine the behaviour of molecules in complex systems. Thus simulation of molecular interaction is essential to predict the results of variable changes that may affect the systems where molecules lie, such as the interactions of a drug with its solvent, or the receptors where it binds. When working with simple molecules standard statistical methods suffice, but in the case of complex molecules such as DNA, proteins or RNA, molecular simulation methods require the inclusion of variables that affect, in apparent serendipity, the global behaviour of biological systems.

Complex simulation can be used, for example, to predict how a protein folds in space, how a membrane receptor transports the ligand to the cellular cytoplasm, how this receptor interacts with the ligand to elicit a physiological response, or how the viral conformation changes to bind to the host membrane and successfully infect it. Complex simulation and prediction systems are in many cases based on the Monte Carlo method, amongst other methods. Furthermore, the Monte Carlo method has been used to simulate and predict evolutionary changes using molecular data from different species.

Further, the folding of chains of protein atoms into shapes that make them usable can be modelled. This modelling can be done using the Monte Carlo method in a simulation (i.e. a Monte Carlo simulation), where each time step simulates a random step forward, and then the equations that create the point of electro-magnetic equilibrium are solved as is known in the art. These equations are extremely complex to calculate and use a significant amount of computing resource.

In structural engineering, for example, for a car manufacturer, there are generally very complex and accurate models that can be used to determine the stress and damage that each component of the car may receive. It is desirable to understand how long a car would last without a component failing in different driving conditions in order to optimise the design and interaction of the components. Different driving scenarios can be simulated over different periods of time (e.g. city driving, on the highway, in the countryside, a combination of city and countryside, at different speeds, etc.). The model can then be used to determine the incremental damage to each component of the vehicle until the damage reaches a threshold of stress that would cause the component to fail. By running these simulations it can be understood how long, on average, the car could remain roadworthy, the probability of it breaking in 1 year, in 5 years, in 20 years, etc.

Typical Prior Art Forecasting System

In the examples described above, the underlying variables that affect the desired parameter may be determined using stochastic modelling such as the Monte Carlo method. The stochastic models can then feed into the determination of the desired result, as the complex relationship between the result and the underlying variable factors is known. However, in order to achieve more accurate results, more combinations of the underlying variables must be run. This becomes computationally very expensive, especially when increasing the number of iterations.

A typical prior art forecasting system 10 is shown in FIG. 1. The forecasting system 10 is suitable for applications wherein the underlying variables that affect the desired parameter may be determined using stochastic modelling. The features of the forecasting system 10 are described in general terms here, and in example contextual applications below. The forecasting system 10 comprises a controller 12, an input factor module 14, a scenario module 16, a metric module 18 and an input/output module 20. The input factor module 14, the scenario module 16, the metric module 18 and the input/output module 20 are each connected to the controller 12.

The input factor module 14 determines the different scenarios 21 (i.e. stochastic models) of the underlying variables that feed into parameter functions 25 (discussed in greater detail below). The input parameter scenarios 21 are stored in a stochastic model database 22.

The metric module 18 is arranged to calculate metrics 19 to generate a forecast based on the evaluated parameter functions 25. The input/output module 20 provides an interface with the forecasting system 10 to external parties. For example, values of the input variables may be provided by a user or an independent system via the input/output module 20.

FIG. 2 shows the scenario module 16 in greater detail. The scenario module 16 comprises a parameter function evaluator 24 and a parameter function database 26 each connected to a scenario controller 28. The parameter function evaluator 24 is arranged to evaluate parameter functions 25 using input variables from the input factor module 14 to generate forecast scenarios 27.

The parameter function database 26 comprises a plurality of application-specific parameter functions 25. Inputs to the parameter function are a set of one or more variables from the input factor module 14. When the parameter function is evaluated for values of the set of one or more variables, the output of the parameter function is a value of the parameter at the given values of those variables.

FIG. 3 shows the relationship between the input scenarios 21, the underlying variables that create the scenarios (in this example, two variables), one of the parameter functions 25 (in this example, function F2), the forecast scenarios 27 and the metrics 19. Namely, each scenario is defined by a different combinations of values of the underlying variables. These scenarios are the input to the function F2 used to generate outputs. The outputs of F2 are the forecast scenarios 27 which in turn can be used to determine metrics 19.

FIG. 4 shows a flowchart summarising the prior art process 30 carried out by the forecasting system 10 described above in relation to FIG. 1. The process 30 begins with all the input variables for the parameter function being evaluated at Step 32 by the input factor module 14 using stochastic modelling and being stored in the stochastic model database 22.

Using the input variables as inputs to the parameter functions 25, the parameter scenarios are determined at Step 34 by the scenario module 16. The parameter scenarios also vary stochastically but are different to the stochastic models of the input variables. The parameter function evaluator 24 of the scenario module 16 evaluates the appropriate parameter function 25 from the parameter function database 26 as many times as required to determine the desired parameter scenarios. The parameter function to be used for the forecast is typically predetermined.

Metrics 19 are then evaluated at Step 36 by the metric module 18 taking input from the parameter scenarios generated by the scenario module 16. The prediction or forecast can then be generated at Step 38 using the metrics determined at Step 36. The output of the forecast may be a probability of an event occurring.

The forecasting system 10 is described in greater detail below in the context of the oil and gas industry, weather forecasting and financial modelling. The oil and gas and weather forecasting examples are described with reference to the typical prior art forecasting system described with reference to FIGS. 1 to 4.

Oil and Gas

In the example of oil and gas industry, the process 30 carried out by the forecasting system 10 involves the input factor module 14 evaluating at Step 32 the stochastic models for input variables. The set of input variables may include: viscosity of the crude oil, remaining supplies in the oil reservoir, oil pumping capacity, well temperature, availability and/or effectiveness of drilling fluid which each vary stochastically over time. The different scenarios 21 (i.e. the stochastic models) for each of the input variables are stored in the stochastic model database 22.

The parameter function database 26 of the scenario module 16 in this example stores a parameter function 25 to determine the extractable volume of crude oil in an oil well over time. The scenarios 21 of this parameter are determined at Step 34 by the scenario module 16 using the input variables from the stochastic model database 22. The parameter function is evaluated at each time step in each scenario to form a stochastic model of the extractable volume parameter.

Metrics are then determined by the metric module 18 at Step 36 based on the expected volume of crude oil that can be extracted from the oil well. In this example, the metrics 19 may include required oil refinement and purification capacity and staffing levels required in downstream processes. The calculated metrics 19 allow the forecast to be generated at Step 38, for example the probability of the future refinement and purification demand exceeding current capacity. Accurate forecasts are beneficial to enable refinement and purification capacity to be increased in a timely manner to meet the predicted demand making the refinement and purification process more efficient.

For the metric of staffing level required, accurate forecasts are advantageous in determining whether additional or fewer staff need to be employed to meet the production demand of the oil well.

In order to improve the accuracy of the forecast, the time resolution of the parameter function scenarios determined at Step 34 is increased. However, the typical time to evaluate the parameter function 25, combined with the number of times the parameter function 25 is evaluated for each scenario 21 means that an accurate high-resolution stochastic model for the extractable volume parameter generally may take a long time to generate even using systems employing the highest power computing devices available. If the time taken to generate the parameter function scenarios is too long, the results may no longer be valid as the value of the starting conditions may have changed.

In order to generate the stochastic model for the extractable volume parameter in a timely manner for the forecast, prior art forecasting systems 10 generally undertake one or more of the following:

1) Reduce the number of scenarios and time points determined at Step 34 by the scenario module 16, so the calculation can be done in a reasonable time. However, the disadvantage is that this reduces the resolution of the scenarios and decreases the accuracy of the metrics used to generate the forecast; 2) Not attempt the calculation of complex parameter functions, as these would utilise more computational time. Again, as a consequence, the determined metrics would be negatively affected by a decrease in accuracy; and 3) Increase the number of computers (and hence computational power) configured to implement the forecasting system 10. However, the number of computers required to reduce the time required to ensure that the results are valid and useful, can be completely impractical. Accordingly, even with increased computing power, the results still suffer from a low resolution and poor accuracy. 4) Implement shortcuts in the methods of the calculation that may increase the speed, but that are always at the expense of accuracy in the calculation. As a result, the accuracy and usability of the metrics generated is compromised.

Weather Forecasting

Another example application of the forecasting system 10 is in the field of weather forecasting. When predicting the temperature or rainfall in a location at a particular time, the temperature or rainfall can be predicted using parameter functions where the set of one or more input variables may include air pressure, air density, solar energy and humidity each from hundreds or more different locations. Monte Carlo simulations are used for weather prediction, starting with the state of the atmosphere when the conditions (e.g. temperature, humidity, pressure) are known.

The process 30 carried out by the forecasting system 10 begins at Step 32, wherein the input factor module 14 determines the variation over time for the input variables across a plurality of scenarios 21. The different scenarios 21 of the variation over time (i.e. the stochastic models) for each of the input variables are stored in the stochastic model database 22.

In this example, the parameter function database 26 includes a parameter function 25 which outputs the temperature at a geographic location as a function of one or more of the input variables. The scenarios 21 of this parameter 25 are determined at Step 34 by the scenario module 16 using the input variables from the stochastic model database 22. The parameter function 25 is evaluated at each time step in each scenario 21 to form a stochastic model of the temperature.

Metrics 19 are then determined by the metric module 18 at Step 36 based on the expected temperature variation at the geographic location. The metrics 19 may include, for example, sea level, pollen count, UV index and wind speed/direction. The calculated metrics 19 allow the forecast to be generated at Step 38. For example, the wind speed and direction metric may be used to determine the probability of the wind speed exceeding a predetermined threshold above which precautionary measures would need to be taken by public services in a timely manner. Another example is if the probability of the pollen count or UV index exceeding a predetermined threshold. The forecast can be publicised to the general public to enable them to take appropriate action.

In weather modelling, as well as in oil and gas, accurate forecasts are advantageous. The accuracy of the metrics 19 is dependent on the resolution of the parameter function scenarios 21 determined at Step 34. However, the typical time to evaluate the parameter function 25, combined with the number of times the parameter function 25 is evaluated for each scenario 21 means that an accurate forecast generally takes a long time to generate. One problem is that if the time taken to generate the parameter function scenarios is too long, the results may no longer be valid as the known value of the starting conditions (from which the input parameters are generated, and that may be included in the input parameters) may have changed.

The parameter function database 26 may also include a parameter function 25 arrange to output the state of equilibrium in the atmosphere, such that a scenario shows the evolution of the atmospheric conditions. Typically, this parameter function is evaluated at Step 34 with small time steps in the Monte Carlo simulation, and after each time step, the state of equilibrium in the atmosphere is determined before the next time step. The smaller the time steps, the easier it is to solve the equations that describe the new state of atmospheric equilibrium. However, with a smaller time step, it becomes computationally expensive to run the simulation far into the future (i.e. a longer time frame). Therefore, to be able to predict further in time in the simulation, larger time steps are used in the simulation. However, these larger time steps make the equations that calculate the new atmospheric equilibrium point more complex and therefore more difficult to solve.

In order to generate the parameter function scenarios in a timely manner for the forecast, prior art forecasting systems 10 generally undertake one or more of the three options expressed above, namely:

1) Reducing the number of scenarios and time points determined at Step 34 by the scenario module 16; 2) Simply not attempting the calculation of complex parameter functions; and 3) Increasing the number of computers (and hence computational power) configured to implement the forecasting system 10. 4) Implementing shortcuts in the methods implemented that may increase the speed, but that decrease the accuracy of the calculation.

Often, none of the above lead to a satisfactory solution that is as accurate as desired.

Financial Modelling

A more detailed but functional example of forecasting is discussed below in the field of financial modelling to illustrate the limitations of existing state of the art forecasting simulation technology. However, it is to be appreciated that the technical limitations which exist in the implementation of these systems and techniques apply equally to all fields of application of forecasting and simulation technology.

Taking finance as a specific example, Over-The-Counter (OTC) derivatives are significant part of the world of global finance. OTC derivatives are contracts signed by two institutions by which they agree to exchange a series of cash-flows during and/or at the end of the contract. Those cash-flows will be driven by the value of some specified external variables. Both institutions agreeing to these contracts are typically called the tounterparties' of the OTC derivative. A financial institution may comprise a large number of OTC derivatives in a ‘portfolio’ or ‘book’, with typically up to of a few million OTC derivatives.

Each OTC derivative comprises one or more price parameters. Examples of price parameters include: a value of the price, a time validity of the price and a potential variation range (e.g. a standard deviation) of the price, a sensitivity of the price to external factors (e.g., the so-called “Greeks” in the art).

An example of a contract is an Equity Future. In this example, today is 1 Jan. 2014 and entity A and entity B agree to the following cash-flow between them on 1 Jan. 2015: the difference in the value of SP500 on those two days, multiplied by an agreed notional value. Accordingly, they agree that a price parameter, P, of the contract is approximately determined by:

P=Notional(SP500_(1jan15) −SP500_(1jan14))

where Notional is, for example, one thousand US Dollars, and SP500 is the value of the Standard & Poor's 500 equity index such that SP5001_(jan14) is the value on 1 Jan. 2014 and SP500_(1jan15) is the value on 1 Jan. 2015. If P is greater than 0, then entity A receives the cash-flow from entity B, if P is less than 0, then entity A pays the cash-flow to entity B.

It must be noted that the price parameter and risk management of these assets, which have a value that changes constantly and, often, very rapidly, is a major issue for asset owners. The reason for this is that the evaluation of accurate risk metrics can be highly computationally demanding to the extent that, in many cases, it cannot be carried out properly. This will be explained in more detail in this document. Furthermore, calculation of price parameters and risk metrics demands the valuation directly from market data using price parameter functions.

This contract has a price parameter at any time, t, between 1 Jan. 2014 and 1 Jan. 2015, which can be approximated by:

P _(t)=Notional(SP500_(t) −SP500_(1jan14))

Generalizing the example above, in order to calculate the price parameter of an OTC derivative at any point in time, there are two distinguishable components:

1) Risk Factors

The price parameter of an OTC derivative is based in the value of several risk factors. In the simplified example provided above, there was only one risk factor, the value of the SP500 equity index. However, in most OTC derivatives, there are a substantially greater number of variable risk factors. These risk factors include several types of interest rates, foreign exchange rates, equity prices, credit spreads, commodity prices, bond prices, etc. A financial institution's book of OTC derivatives will be sensitive to a plurality of these factors, typically ranging from a few hundred to several thousand risk factors.

2) Derivative Price Parameter Functions

Once the risk factors that affect the price parameter of the OTC derivative have been determined, the price parameter at time t will be given by a function ƒ(x), such that:

P _(t)=ƒ(x _(t))

where x refers to the weighted and combined collection of all the risk factors that are relevant to the price parameter of that OTC derivative at time t. Referring back to the previous example, if x is taken to be equivalent to SP500, then ƒ(x_(t))=Notional (x_(t)−x₀), where x₀ is the initial value of the risk factor, i.e. SP500_(1jan14).

Determination of Counterparty Credit Risk Metrics

In the context of financial services, counterparty credit risk is the risk that a counterparty to a financial contract will default prior to the expiration of the contract. As OTC derivatives are generally privately-negotiated contracts between counterparties, they are subject to counterparty risk. Counterparty risk is similar to other forms of credit risk in that the cause of economic loss is if the obligor defaults (e.g. goes into liquidation). However, counterparty risk is differentiated from more conventional forms of credit risk in the uncertainty of exposure (i.e. the value of the potential loss may vary) and its bilateral nature (i.e. both parties are exposed to each other). In order to reduce the uncertainty of exposure and manage the counterparty credit risk, the risk is measured and quantified using one or more metrics.

FIG. 5A shows an example prior art forecasting system 100 that can be used to determine counterparty credit risk metrics. The system comprises a controller 102 and a risk factor module 104, a price parameter module 106, a risk metric module 108 and an input/output module 110 that are each connected to the controller 102. The risk factor module 104 determines the variation over time of the underlying variables that feed into the price parameter functions and affect the price of the OTC derivatives.

A price parameter function, ƒ(x), is a mathematical function that defines a relationship between one or more input variables to produce an output, where the output is a parameter of the price. Price parameter functions are defined in the domain(s) of the one or more input variables. Accordingly, a price parameter function varies in one or more relationship domains (as opposed to the time domain) where each relationship domain corresponds to the effect of the one or more input variables on the output of the function.

The risk factor module 104 stores data regarding the value of the risk factors in a stochastic model database 112, described in more detail with reference to FIG. 5B. The price parameter module 106 determines the price parameter function ƒ(x), and is described in more detail with reference to FIG. 6A. The risk metric module 108 is configured to determine the value of the counterparty credit risk metrics. The input/output module 110 provides an interface with the forecasting system 100 to external parties.

Each module comprises its own processing core or set of processors arranged in parallel with a multicore structure to provide increased processing speed as is known in the art. The input/output module 110 is configured to allow two-way communication between the forecasting system 100 and a user (user device or computer) or computer network, wherein inputs are signals or data received by the system, and outputs are signals or data sent from it. For example, the input/output module 110 may be a network card, a modem or a computer terminal comprising a keyboard and a monitor.

FIG. 5B shows the stochastic model database 112 in greater detail. For each risk factor 114, for a point in time 116, values 120 of the risk factor are stored for a plurality of different scenarios 118. In the example of FIG. 5B, the values 120 of two different risk factors, Factor1 and Factor2, are shown each at four time points.

FIG. 5C is a graphical representation of an example stochastic model 122 that may be stored in the stochastic model database 112. The stochastic model 122 shows the value of a risk factor over time with four distinct distributions of potential outcomes 124, called ‘scenarios’. In the example shown all four scenarios 124 start at t=0 with the same value of the risk factor; however it must be noted that, in other applications, this may not always the case and the initial value of the scenarios may be different. In FIG. 5C, only four scenarios are shown for illustrative simplicity, however, a typical stochastic model may comprise thousands of scenarios.

In this example, the initial conditions of each risk factor are known, in other words, the value of each risk factor is known for the present time, t=0, however the values of the risk factor any time after the present time are not known. As each risk factor may vary over time, the risk factor module is configured to determine the evolution of the risk factors over time using probabilistic modelling in the form of stochastic simulation such as the Monte Carlo method as is known in the art.

The four scenarios 124 tend to diverge away from each other with increasing time as the stochastic model allows for random variations which accumulate, causing increased uncertainty in the potential future value of the risk factor. The random variations can be based on fluctuations observed in historical data of the risk factor and/or in the value of current data.

As described above, each OTC derivative may be affected by a few hundred to several thousand risk factors, each with their own stochastic model of value for thousands of scenarios.

The stochastic scenarios generated by the risk factor module 104 are communicated to the price parameter module 106 via the controller 102.

The price parameter module 106 is configured to determine the price parameter of each OTC derivative at a desired time point and in a given scenario. For a particular time, t, and in a given scenario, i, the price parameter of the OTC derivative can be determined using a pricing parameter function, P_(t,i)=ƒ (x_(t,i)). To establish the variation of the OTC derivative price parameter over time in a given scenario, the pricing parameter function can be determined for different points in time.

FIG. 6A shows an example prior art price parameter module 106 in greater detail. The price parameter module 106 comprises a price parameter controller 150, a price parameter function evaluator 152 and a pricing function database 154. The price parameter controller 150 is connected to the controller 102 of the forecasting system 100. The price parameter function evaluator 152 and the pricing function database 154 are each connected to the price parameter controller 150. The price parameter function evaluator 152 determines the value of a price parameter function from the pricing function database 154 using the value of the risk factors from the stochastic model database 112. The pricing function database 154 comprises a plurality of price parameter functions.

FIG. 6B shows the pricing function database 154 in greater detail. The pricing function database 154 stores the price parameter function 156 associated with each OTC derivative 158 in the financial institution's portfolio. Each price parameter function 156 generally comprises a different combination and weighting of risk factors. The price parameter function evaluator 152 is configured to determine the output value of the price parameter function 156, requesting the value of the risk factors it requires from the stochastic model database 112 via the risk factor module 104 and the controller of the system and the price parameter controller 150.

FIG. 7A is a graphical representation 180 of an example of a price parameter function over time and shows an example distribution in the price parameter of an OTC derivative over a period of time in two potential outcomes 182 (i.e. two scenarios). The price parameters of the OTC derivative in each of the two scenarios is determined by the price parameter function evaluator 152 and are shown at each of ten time points, t₀ to t₉ respectively.

FIG. 7B is a graphical representation 184 of another example of a price parameter function over time showing a higher resolution distribution of four scenarios 186 compared to the lower resolution shown in FIG. 7A. The higher resolution means that the price parameter function has been evaluated at more time points.

Increasing the number of time points at which the price parameter is determined increases the resolution of the scenario, thereby increasing the accuracy of any analysis made on the scenario. However, as is described below, this comes at a far greater computational cost. This is because each determination of the price parameter function by the price parameter module 106 is computationally expensive, and the price parameter function must be evaluated at each desired time point and scenario.

The number of scenarios in the stochastic model generated by the risk factor module 104 will be referred to below as N. For example, referring to FIG. 5C, N=4. Additionally, the number of points in time at which the OTC price parameter is determined at will be referred to M. For example, referring to FIG. 7A, M=M. It follows that in a price parameter distribution, the price parameter function is determined N×M times for each derivative. For example, referring to FIG. 7A, the price parameter function is called by the price parameter function evaluator 2×10=20 times.

The price parameter distributions generated by the price parameter module 106 are communicated to the risk metric module 108 via the controller 102.

The risk metric module 108 is configured to determine a plurality of counterparty credit risk metrics using all the scenarios in the price parameter distributions for all the OTC derivatives of the financial institution. Many metrics are typically used for pricing, accounting, risk management and capital calculations. The metrics include:

a) Expected Exposure (EE) profiles, both in its ‘positive’ (EPE) and ‘negative’ (ENE) version b) Potential Future Exposure (PFE) profiles c) Credit Value Adjustment (CVA) and Debit Value Adjustment (DVA) for pricing d) Effective Expected Positive Exposure (EEPE) for regulatory capital calculation e) Funding Value Adjustment (FVA) for internal valuation and balance sheet calculation.

FIG. 8A shows a flowchart summarising the prior art process 200 carried out by the forecasting system 100 described above in relation to FIG. 5A. The process 200 begins at Step 202, where all the risk factors are evaluated using stochastic modelling by the risk factor module 104 and stored in the stochastic model database 112. Using risk factors as inputs to the price parameter functions, the price parameter scenarios (which are related but different to the risk factor scenarios) are determined at Step 204 by the price parameter module 106. The price parameter function evaluator 152 of the price parameter module 106 evaluates the appropriate price parameter function from the pricing function database 154 as many times as required to determine the desired price parameter scenarios (described in more detail with reference to FIG. 8B). The risk metrics are then evaluated at Step 206 by the risk metric module 108 taking input from the price parameter scenarios generated by the price parameter module 106. The prediction or forecast can then be generated at Step 208 using the risk metrics determined at Step 206.

FIG. 8B shows a flowchart of the process 204 carried out by the price parameter module 106 in order to determine the price parameter scenarios. The process 204 starts with the price parameter function evaluator 152 receiving at Step 210 the number of time points, M, required for the desired resolution of the price parameter scenarios. Then the price parameter function evaluator 152 receives at Step 212 the number of scenarios, N, desired. The price parameter function evaluator 152 then obtains at Step 214 the price parameter function 156, ƒ(x), for the OTC derivative from the pricing function database 154.

The price parameter function evaluator 152 sets at Step 216 the current time point to the first time point of the total number of time points, M. Then, the price parameter function evaluator 152 sets at Step 218 the current scenario to the first scenario of the total number of scenarios, N.

For the first time point of the first scenario, the price parameter function 156 is then evaluated at Step 220 by the price parameter function evaluator 152. The price parameter function evaluator 152 then checks at Step 222 whether the price parameter function has been evaluated for all scenarios N received at Step 212 for the first time point. If not, the price parameter function evaluator 152 moves at Step 224 the current scenario to the next scenario and the process 204 loops back to Step 220 whereby the price parameter function 156 is evaluated by the price parameter function evaluator 152 for the current time point.

This loop continues until the price parameter function 156 has been evaluated for all the scenarios, N, received at Step 212 for the first time point. If, following the check at Step 222, the price parameter function 156 has been evaluated for all the scenarios received at Step 212, the process 204 moves on to checking at Step 226 whether price parameter function 156 has been evaluated for all time points received at Step 210.

If not, the price parameter function evaluator 152 moves at Step 228 the current time point to the next time point and then the current scenario is moved at Step 230 to the first scenario. The process 204 then loops back to Step 220 whereby the price parameter function 156 is evaluated by the price parameter function evaluator 152 for the current time point of the current scenario. Once the price parameter function 156 has been evaluated for all the time points for all the scenarios, the process 204 is completed.

In other embodiments, all the time points for the first scenario are determined before moving on to determining all the time points for the next scenario, until all the time points for all the scenarios have been determined.

The process 204 carried out by the price parameter module 106 in order to determine the price scenarios demonstrates that the number of evaluations of the price parameter function 156 required (M×N in this example) is a limiting factor in the speed of the overall process 200 carried out by the forecasting system 100. A more detailed numerical example of the limitations imposed by the process 204 is highlighted in the financial example below.

Example Number of Price Parameter Function Determinations

Taking the CVA metric as an example, the price parameter of the OTC derivative is typically desired to be calculated to an accuracy of one ‘basis point’ (one one-hundredth of one percent), i.e. 0.01% ( 1/10⁴). The typical number of scenarios, N, used in order to achieve this pricing error in CVA is around one million.

In order to calculate the CVA ‘Greeks’ via Monte Carlo methods (the sensitivities of the CVA price parameter to changes in the market conditions), the CVA price parameters must be calculated to a resolution of 0.01 basis points at least ( 1/10⁶) so that the calculated Greek parameter has a stable value. As the pricing error expands in Monte Carlo simulations as √{square root over (N)}, this means that in order to calculate the Greeks to this accuracy, N would be around ten billion (i.e. 10,000,000,000) and this would be received by the price parameter function evaluator 152 at Step 212 of the process 204.

The expiry of an OTC derivative may range from a few weeks (e.g. foreign exchange derivatives) to up to 50 years (e.g. inflation derivatives). The average expiry in a financial institution's book of derivatives can typically be of around 10 years.

For example, taking an average OTC derivative with an expiry of 10 years, there are approximately 2,600 trading days in that period. Ideally, the financial institution would desire an accurate determination of the counterparty credit risk that it is currently exposed to be based on the forecast from today (i.e. to) to the end of year 10. Taking each day as a point in time at which the price parameter is determined for one of the OTC derivatives, M=2,600 and this would be received by the price parameter function evaluator 152 at Step 210 of the process 204. Therefore, in order calculate, for example, the CVA counterparty credit risk metric for this OTC derivative to the desirable degree of accuracy of 0.01 basis points, the price parameter function is determined N×M times or 10,000,000,000×2,600=26,000,000,000,000 (i.e. 26×10¹² times) in the step of determining the price parameter scenarios (Step 220 of the process 204 carried out by the price parameter function evaluator 152).

A typical large financial institution may have a few million OTC derivatives on their books. For example, if the number of derivatives, D, of a financial institution is 1,000,000, then to calculate the CVA metric for each (and assuming an average expiry of 10 years for each OTC derivative), the number of times that the price parameter function 156 is determined is D×N×M times or 1,000,000×26,000,000,000,000=26,000,000,000,000,000,000 (i.e. 26×10¹⁸ times).

This number of calculations would only determine one set of risk metrics. Additionally, a number of subsequent calculations are carried out. In order to calculate the CVA Greeks via the typical Monte Carlo method, for example, this calculation would be repeated at least one time per Greek. A large financial institution may typically desire to calculate a few hundred Greeks. The example below uses 100 Greeks.

Also, whenever a new derivative is transacted, an existing derivative is unwound, or when risk analysis is performed, the calculation would be to be re-run, at least partially to the netting set affected by the derivative.

Conservatively, a typical financial institution may desirably run this calculation 10 times per day.

This would lead to the number of price parameter function determinations, in this example, to be:

100×10×26×10¹⁸=26×10²¹

The number of calculations required to achieve this desired level of accuracy demands significant computational effort and accordingly, takes too long to process to be of any benefit. Accordingly, this level of accuracy is simply not achievable within a practical time period given the available computing power.

Maximising the accuracy of the risk metrics is desirable in all fields, including weather forecasting and in the oil and gas industry. In turn, this increases the number of times the respective functions ƒ(x) are evaluated, requiring more computational effort and time.

Computational Power Requirements

The largest financial institutions in the world comprise server farms of state-of-the-art computers, available 24 hours a day, for these calculations. In spite of that, the reported number of scenarios that they use range between N=1,000 and N=10,000, and the number of time points between M=25 to M=100.

The computational time required for an optimised price parameter function 156 for a single parameter currently range from around one millisecond for the simplest derivatives to several minutes or hours for complex ones. Reportedly, large financial institutions typically run a few hundred systems, typically about two hundred, constantly and simultaneously to determine the price parameter functions and process the counterparty credit risk metrics.

The approximate total computational time of Step 204 is given by:

computational time per calculation×number of derivatives,D×number of scenarios,N×number of time points,M number of computers

Assuming that an average computational time for a price parameter function is only 0.01 seconds, taking the typical values above, this gives:

$\mspace{20mu} {\frac{0.01 \times 1,000,000 \times 1,000 \times 25}{200} = {{1,250,000\mspace{14mu} {seconds}} \approx {14.5\mspace{14mu} {days}}}}$ $\frac{0.01 \times 1,000,000 \times 10,000 \times 100}{200} = {{50,000,000\mspace{14mu} {seconds}} \approx {578.7\mspace{14mu} {days}}}$

Therefore the computational time of Step 204 can easily range from 14.5 days to 578.7 days.

The vast number of determinations of the price parameter function 156 by the price parameter function evaluator 152 means that by the time the calculations are complete, the results may no longer be valid or useful. Therefore, financial institutions implement the forecasting system 100 in state-of-the-art computer farms but generally reduce the accuracy in order to save time. Further, in order to achieve the best possible result in the available time, workarounds are implemented including:

1) Reducing the number of scenarios and time points so the calculation can be done in a reasonable time. As a result, the quality of the risk metrics and subsequent risk management results generated by the forecasting system 100 substantially decreases. 2) Simply not attempting this calculation for sophisticated OTC derivatives with complex price parameter functions, as these would utilise more computational time. Again, as a consequence, risk management results are negatively affected. 3) Increasing the number of computers (and hence computational power) configured to implement the forecasting system 100. However, the number of computers required to reduce the time required to ensure that the results are valid and useful is impractical. Further, the size of the computational problem is so large that additional computing power can actually make a negligible difference to the time required. 4) Implement methodological shortcuts in the calculation that may increase the speed, but that always create lack of accuracy in the calculation and, hence, decrease the usability of the metrics produced.

The present invention seeks to mitigate or at least overcome at least some of the above-mentioned problems.

SUMMARY OF THE INVENTION

According to the present invention, there is provided a computer-implemented method of determining an approximated value of a parameter in a first domain, the parameter being dependent on one or more variables which vary in a second domain, and the parameter being determined by a function which relates sets of values of the one or more variables in the second domain to corresponding values in the first domain, the method being implemented on a computer system including a processor, and the method comprising: determining a plurality of anchor points in the second domain, wherein each anchor point comprises a set of values of the one or more variables in the second domain; evaluating, at each anchor point, the function to generate corresponding values of the parameter in the first domain; generating an approximation function to the function by fitting a series of orthogonal functions or an approximation to a series of orthogonal functions to the corresponding values of the parameter in the first domain; and using the approximation function to generate the approximated value of the parameter in the first domain.

Optionally, the method comprises transforming the function into an intermediary function before the generating step, wherein the intermediary function is a linear transformation of the function, and wherein the generating step comprises approximating the intermediary function.

Optionally, the transforming step comprises linearly transforming the function into a predetermined range.

Optionally, the series of orthogonal functions is one of the group comprising: sines, cosines, sines and cosines, Bessel functions, Gegenbauer polynomials, Hermite polynomials, Laguerre polynomials, Chebyshev polynomials, Jacobi polynomials, Spherical harmonics, Walsh functions, Legendre polynomials, Zernike polynomials, Wilson polynomials, Meixner-Pollaczek polynomials, continuous Hahn polynomials, continuous dual Hahn polynomials, a classical polynomials described by the Askey scheme, Askey-Wilson polynomials, Racah polynomials, dual Hahn polynomials, Meixner polynomials, piecewise constant interpolants, linear interpolants, polynomial interpolants, gaussian process based interpolants, spline interpolants, barycentric interpolants, Krawtchouk polynomials, Charlier polynomials, sieved ultraspherical polynomials, sieved Jacobi polynomials, sieved Pollaczek polynomials, Rational interpolants, Trigonometric interpolants, Hermite interpolants, Cubic interpolants, and Rogers-Szegö polynomials.

Optionally, the series of orthogonal functions is a series of orthogonal polynomials.

Optionally, the approximation to a series of orthogonal functions is an approximation to a series of orthogonal polynomials.

Optionally, the method comprises selecting a series of orthogonal functions or an approximation to a series of orthogonal functions before the generating step, and wherein the generating step comprises using the selected series of orthogonal functions or the approximation to the series of orthogonal functions.

Optionally, the number of anchor points, N_(E), is so that:

$\frac{N_{S} \cdot T}{{N_{E} \cdot T} + {N_{S} \cdot t}} > 1$

Wherein N_(S) is the number of scenarios, t is the time taken to run the approximation function, and T is the time taken to run the function. For example, the number of anchor points is from 1 to 200, from 5 to 50, from 1 to 50, or from 5 to 20.

Optionally, the method comprises holding the values of all but one variable in the set to be constant in the evaluating step.

Optionally, the method comprises using the approximated values of the parameter to generate standard metrics.

Optionally, the values of the variables vary stochastically.

Optionally, the using step comprises using the approximation function a plurality of times to generate at least one scenario or one time step, each scenario or each time step comprising a plurality of approximated values.

Optionally, the number of scenarios or time steps is significantly greater than the number of anchor points.

Optionally, an output of the function is a parameter of a financial product.

Optionally, the financial product is a financial derivative including one or more of a group comprising: an option pricing function, a swap pricing function and a combination thereof.

Optionally, the function is one of the group comprising: a Black-Scholes model, a Longstaff-Schwartz model, a binomial options pricing model, a Black model, a Garman-Kohlhagen model, a Vanna-Volga method, a Chen model, a Merton's model, a Vasicek model, a Rendleman-Bartter model, a Cox-Ingersoll-Ross model, a Ho-Lee model, a Hull-White model, a Black-Derman-Toy model, a Black-Karasinski model, a Heston model, a Monte Carlo based pricing model, a binomial pricing model, a trinomial pricing model, a tree based pricing model, a finite-difference based pricing model, a Heath-Jarrow-Morton model, a variance gamma model, a Fuzzy pay-off method, a Single-index model, a Chepakovich valuation model, a Markov switching multifractal, a Datar Mathews method, and a Kalotay-Williams-Fabozzi model.

Optionally, the anchor points are the zeros of each of the orthogonal functions, or a subset of the zeros.

Optionally, the approximating function is generated using an interpolation scheme.

Optionally, the interpolation scheme is one from the group comprising: piecewise constant interpolants, linear interpolants, polynomial interpolants, gaussian process based interpolants, spline interpolants, barycentric interpolants, Rational interpolants, Trigonometric interpolants, Hermite interpolants and Cubic interpolants.

Optionally, the anchor points are integration points of a numerical integration scheme.

Optionally, the numerical integration scheme is one from the group comprising: Newton-Cotes methods, a trapezoidal method, a Simpson's method, a Boole's method, a Romberg integration method, Gaussian quadrature methods, Chenshaw-Curtis methods, a Fejer method, a Gaus-Kronrod method, Fourier Transform methods, an Adaptive quadrature method, a Richardson extrapolation, a Monte Carlo and Quasi Monte Carlo method, a Markov chain Monte Carlo, a Metropolis Hastings algorithm, a Gibbs Sampling, and Fast Fourier Transform methods.

Optionally, the speed of the calculation is increased with no loss of accuracy in the standard metrics, and/or the accuracy of the standard metrics increase with no decrease in the speed of the calculation when using the approximation function.

Optionally, there is a difference between the approximated value of the parameter and the parameter in the 2^(nd) significant figure when compared at the same point in the first domain.

Optionally, there is a difference between the approximated value of the parameter and the parameter in between the 4^(th) and the 6^(th) significant figure when compared at the same point in the first domain.

Optionally, there is a difference between the approximated value of the parameter and the parameter in the 15^(th) significant figure when compared at the same point in the first domain.

A financial derivative may comprise a parameter, wherein a value of the parameter is determined using the computer-implemented method above.

According to an aspect of the invention, there is provided a computer system comprising a processor configured to determine an approximated value of a parameter in a first domain, the parameter being dependent on one or more variables which vary in a second domain, and the parameter being determined by a function which relates sets of values of the one or more variables in the second domain to corresponding values in the first domain, the computer system comprising: a determination module arranged to determine a plurality of anchor points in the second domain, wherein each anchor point comprises a set of values of the one or more variables in the second domain; an evaluation module arranged to evaluate, at each anchor point, the function to generate corresponding values of the parameter in the first domain; a generation module arranged to generate an approximation function to the function by fitting a series of orthogonal functions or an approximation to a series of orthogonal functions to the corresponding values of the parameter in the first domain; and an approximation module arranged to use the approximation function to generate the approximated value of the parameter in the first domain.

According to an aspect of the invention, there is provided a computer-implemented method of determining an approximated value of a parameter in a first domain, the parameter being dependent on one or more variables which vary in a second domain, and the parameter being determined by a function which relates sets of values of the one or more variables in the second domain to corresponding values in the first domain, the method using a series of orthogonal functions or an approximation to a series of orthogonal functions to approximate the function.

BRIEF DESCRIPTION OF THE DRAWINGS

FIGS. 1 through 8B have already been described above by way of background, in which:

FIG. 1 is a schematic block diagram of a prior art forecasting system;

FIG. 2 is a schematic block diagram of a scenario module of the system of FIG. 1;

FIG. 3 is a schematic diagram of data flow in a prior art forecasting system;

FIG. 4 is a flowchart showing a process carried out by the prior art system of FIG. 1;

FIG. 5A is a schematic block diagram of a prior art forecasting system;

FIG. 5B is a schematic block diagram of a stochastic model database of the prior art system of FIG. 5A;

FIG. 5C is an example chart showing the variation of values of risk factor over time in a stochastic model stored in the stochastic model database of FIG. 5B;

FIG. 6A is a schematic block diagram of a price parameter module of the system of FIG. 5A;

FIG. 6B is a schematic block diagram of a pricing function database of the price parameter module of FIG. 6A;

FIGS. 7A and 7B are example charts showing the variation of price parameters determined at specific time points;

FIG. 8A is a flowchart showing a process carried out by the prior art system of FIG. 1A;

FIG. 8B is a flowchart showing a process of determining price scenarios in the process of FIG. 8A;

One or more embodiments of the invention will now be described, by way of example only, with reference to the remaining drawings, in which:

FIG. 9A is a schematic block diagram of a forecasting system according to an embodiment of the invention;

FIG. 9B is a schematic block diagram of an approximation module of the forecasting system of FIG. 9A;

FIG. 10A is a flowchart showing a process carried out by the forecasting system of FIG. 9A according to an embodiment of the invention;

FIG. 10B is a flowchart showing a process of determining price scenarios in the process of FIG. 10A;

FIG. 10C is a flowchart showing a process of generating approximating functions in the process of FIG. 10B;

FIGS. 11A and 11B are graphical charts comparing obtained Expected Positive Exposure profiles of the prior art forecasting system to the embodiment of FIG. 9A;

FIGS. 12A and 12B are graphical charts comparing obtained Potential Future Exposure profiles of the prior art forecasting system to the embodiment of FIG. 9A;

FIG. 13 is a graphical chart comparing computing time profiles of the prior art forecasting system to the embodiment of FIG. 9A;

FIGS. 14A and 14B are graphical charts showing the reduced noise of the results obtained using the embodiment of FIG. 9A compared with the prior art forecasting system; and

FIG. 15 is a graphical chart illustrating an application of the embodiment in more than one dimension.

DETAILED DESCRIPTION

The present invention can be realised in many different ways and as has been explained previously for many different applications. In the following description, a specific embodiment is described which is in the financial forecasting domain. The reason for illustrating the present invention using this embodiment is because the present invention can be described efficiently using some of the elements of a system for forecasting future price of an OTC contract which have already been described in the background section. Also benchmarking results using the prior art system 100 and the forecasting system of the present embodiment can be carried and these results have been provided towards the end of this specific description. However it is to be clearly appreciated that the present embodiment is exemplary only and the structure of the specific embodiment described herein can readily be adapted by the skilled person to work with any of the other areas of application mentioned in the introduction.

As described above, the step of determining the price parameter scenarios in the process 204 described in relation to FIG. 8B comprises evaluating the price parameter function 156 potentially quintillions of times and is the most computationally expensive part of the overall process 200 of determining counterparty credit risk metrics.

FIG. 9A shows a forecasting system 300 according to an embodiment of the present invention. The forecasting system 300 has substantially the same structure as the prior art forecasting system 100 shown in FIG. 1. However, the prior art price parameter module 106 described with reference to FIG. 6A is replaced in the current embodiment with an approximation module 302 to significantly reduce the computational time required to generate the required results of the forecast. The modules and database shown with dashed lines in FIG. 9A are substantially unchanged from the prior art forecasting system 100.

The approximation module 302 is configured to evaluate an approximation to the pricing parameter functions 156, for example by substituting the price parameter function 156 by an expansion in orthogonal functions that yields a substantially equivalent result to that of the original pricing parameter function. The orthogonal functions are generally faster to evaluate than their corresponding original price parameter functions.

FIG. 9B shows the approximation module 302 in greater detail. The approximation module 302 comprises an approximation controller 310, a transformer 312, a pricing function database 154 (substantially identical to the prior art pricing function database shown in FIGS. 2A and 2B), an approximation generator 314, an expansion series database 320, an approximation evaluator 316 and an approximating pricing function database 318. The transformer 312, the approximation generator 314, the approximation evaluator 316, the pricing function database 154 and the approximating pricing functions database 318 are each connected to the approximation controller 310.

The approximation generator 314 is arranged to determine as approximation of the pricing functions 156. The expansion series database 320 is connected to the approximation generator 314. The expansion series database 320 comprises expansion series 321 associated with different families of orthogonal functions, including orthogonal polynomials and approximations to them. For example, the expansion series database 320 stores mathematical relationships for orthogonal functions including sines and cosines, Bessel functions, Gegenbauer polynomials, Hermite polynomials, Laguerre polynomials, Chebyshev polynomials, Jacobi polynomials, Spherical harmonics, Walsh functions, Legendre polynomials, Zernike polynomials, Wilson polynomials, Meixner-Pollaczek polynomials, continuous Hahn polynomials, continuous dual Hahn polynomials, the classical polynomials described by the Askey scheme, Askey-Wilson polynomials, Racah polynomials, dual Hahn polynomials, Meixner polynomials, Krawtchouk polynomials, piecewise constant interpolants, linear interpolants, polynomial interpolants, gaussian process based interpolants, spline interpolants, barycentric interpolants, Charlier polynomials, sieved ultraspherical polynomials, sieved Jacobi polynomials, sieved Pollaczek polynomials, Rational interpolants, Trigonometric interpolants, Hermite interpolants, Cubic interpolants, and Rogers-Szego polynomials which are examples of families of mathematical relationships for functions that are suitable for approximating price parameter functions 156. The orthogonal functions stored in the expansion series database 320 are termed expansion series 321. The expansion series 321 may be orthogonal functions or approximations to orthogonal functions.

Accordingly, the transformer 312 is arranged to transform the price parameter functions 156 into the range where an approximating expansion series is substantially equivalent to the price parameter function, if required, prior to the approximation generator 314 determining an approximation of a pricing function 156.

The approximating pricing function database 318 stores the approximations 319 (i.e. approximating functions) of the pricing functions 156 determined by the approximation generator 314.

The approximation evaluator 316 is configured to determine the approximate value of the price parameter using the approximation 319 (also termed as ‘approximation functions’). The function of the approximation evaluator 316 is roughly analogous to the pricing parameter function evaluator 152 of the prior art price parameter module 106.

FIG. 10A shows an overview of the process 348 carried out by the forecasting system 300. The process begins at Step 202, where all the risk factors are evaluated using stochastic modelling by the risk factor module 104 and stored in the stochastic model database 112 in substantially the same way as in the prior art forecasting system 100. Then, the approximation evaluator 316 evaluates at Step 354 the price scenarios using the approximating function(s) 319 from the approximating pricing function database 318. The specific embodiment is described below with one approximating function 319 generated for each time step (Step 354 is described in more detail with reference to FIG. 10B), however, it is to be appreciated that other embodiments are possible which may comprise more than one approximating function 319 generated for each time step.

The risk metrics are then evaluated at Step 206 by the risk metric module 108 taking input from the price scenarios generated by the approximation module 302 in substantially the same way as in the prior art forecasting system 100. The prediction or forecast can then be generated at Step 208 using the risk metrics determined at Step 206 in substantially the same way as in the prior art forecasting system 100.

FIG. 10B shows the process of determining at Step 354 the price scenarios using the approximating function 319 carried out by the approximation module 302 in greater detail.

The process 354 starts with the approximation evaluator 316 receiving at Step 210 the number of time points, M, required for the desired resolution of the scenarios in substantially the same way as in the prior art forecasting system 100. Then the approximation evaluator 316 receives at Step 212 the number of scenarios, N, desired in substantially the same way as in the prior art forecasting system 100. The approximation evaluator 316 then obtains at Step 214 the price parameter function 156, ƒ(x), for the OTC derivative from the pricing function database 154 desired in substantially the same way as in the prior art forecasting system 100.

The approximation evaluator 316 sets at Step 216 the current time point to the first time point of the total number of time points, M in substantially the same way as in the prior art forecasting system 100. Then, the approximation evaluator 316 sets at Step 218 the current scenario to the first scenario of the total number of scenarios, N in substantially the same way as in the prior art forecasting system 100. Following this, a function approximating the pricing function (i.e. the approximating function 319) for the current time point is generated at Step 380 by the approximation module 302 and stored in the approximating pricing function database 318 (this step is described in more detail with reference to FIG. 10C).

For the first time point of the first scenario, the approximating function 319 is then evaluated at Step 382 by the approximation evaluator 316. The approximation evaluator 316 then checks at Step 222 whether the price parameter function has been evaluated for all scenarios N received at Step 212 for the first time point in substantially the same way as in the prior art forecasting system 100. If not, the approximation evaluator 316 moves at Step 224 the current scenario to the next scenario and the process 354 loops back to Step 382 whereby the approximating function 319 is evaluated by the approximation evaluator 316 for the current time point.

This loop continues until the approximating function 319 has been evaluated for all the scenarios, N, received at Step 212 for the first time point in substantially the same way as in the prior art forecasting system 100. If, following the check at Step 222, the approximating function 319 has been evaluated for all the scenarios received at Step 212, the process 354 moves on to checking at Step 226 whether approximating function 319 has been evaluated for all time points received at Step 210 in substantially the same way as in the prior art forecasting system 100.

If not, the approximation evaluator 316 moves at Step 228 the current time point to the next time point and then the current scenario is moved at Step 230 to the first scenario in substantially the same way as in the prior art forecasting system 100. The process 354 then loops back to Step 380 whereby a function approximating the price parameter function is generated for the current time point. Then, at Step 382 the approximating function 319 is evaluated by the approximation evaluator 316 for the current time point of the current scenario. Once the approximating functions 319 have been evaluated for all the time points for all the scenarios, the process 354 is completed.

In other embodiments, all the approximating functions 319 for all the time points, M, are generated at the start of the process. Then the price parameters for all time points for the first scenario are determined before moving on to determining them for all the time points for the next scenario, until all the time points for all the scenarios have been determined. Also, in other embodiments, the loop in time steps and scenarios could be inverted as it was explained in the prior art system 100.

As can be seen, the main difference between the Step 354 carried out the forecasting system 300 over the Step 204 carried out by the prior art forecasting system 100 is that the approximation evaluator 316 evaluates a approximating function 319 which is computationally cheap to evaluate compared to the price parameter function evaluator 152 which evaluates price parameter functions 156 directly at a far higher computational cost.

FIG. 10C shows the process of generating at Step 380 the approximating function 319 carried out by the approximation module 302 in greater detail. The process starts with all but one of the factors in the pricing function being fixed at Step 360 by approximation generator 314 such that all but one degree of freedom (DOE) of the pricing function are locked.

The approximation controller 310 then determines at Step 362 whether the variable, x, in which the price parameter function range is defined over the same range as the approximation function. If not, then the transformer 312 carries out at Step 364 a linear transformation to the variable x of the price parameter function 156 from the pricing function database 154.

If no transformation is required following the check of Step 362 or once the pricing function has been transformed at Step 364, the approximation generator 314 determines at Step 366 a function that approximates the linearly transformed price parameter function using the expansion series 321 from the expansion series database 320. The approximation generator 314 stores at Step 368 the resulting approximation function 319 that approximates the linearly transformed price parameter function (i.e. orthogonal approximations 319) in the approximating pricing function database 318. As will be discussed in more detail below, the polynomial approximations 319 yield substantially the same result as the original price parameter functions 156, but in a significantly reduced time frame.

The steps in the process 348 carried out by the approximation module 302 are described in more detail in the specific example given below, where Chebyshev polynomials are employed by the approximation module 302. However, it is to be appreciated that any family of orthogonal functions or approximations to them can be used to approximate the price parameter functions 156, each with a slightly different degree of accuracy and amounts of computational cost. However, regardless of which function 319 is used, the reduction in computing overhead compared to the prior art forecasting system 100 can, in most cases, be orders of magnitude less.

In other embodiments, the time variable in the price parameter function can be considered as another variable x in a multidimensional generation of the approximation function 319 (explained in greater detail later), in which case it is not necessary to generate a new approximation function for each time step in 380.

The transformer 312 is configured to perform a linear transformation of the price parameter function from the pricing function database, ƒ(x) in which x varies between a minimum value, a, and a maximum value, b. Accordingly, the transformer 312 determines a linearly transformed price parameter function g(y) in which y varies between a minimum value, c, and a maximum value, d. In mathematical notation, this would be:

$\begin{matrix} {x \in \left. \left\lbrack {a,b} \right\rbrack\rightarrow y \right. \in \left\lbrack {c,d} \right\rbrack} \\ \left. {f(x)}\rightarrow{g(y)} \right. \end{matrix}$

In most cases, the transformer 312 carries out Step 364 on the price parameter function and on the variable x for which the price parameter function is evaluated at Step 380. The interval [a,b] is driven by the interval of the original price parameter function in the x domain in which accurate approximation is desired. Typically, this interval is defined by the minimum and the maximum values of the variable x.

The interval [c,d] is defined by the family of orthogonal functions chosen for the approximation. The interval is driven by the interval in the y domain where the approximation is accurate. For example, in the case of Chebyshev polynomials, the interval is [−1,1], and in the case of {sin(nx) and cos(mx)} orthogonal functions, the interval is [−π,π].

Once the transformer 312 has determined the linearly transformed pricing parameter function, the approximation generator 314 can then generate at Step 366 a function that approximates the linearly transformed pricing parameter function.

Using the expansion series 321 in the expansion series database 320 and a price parameter function 156 from the pricing function database 154, the approximation generator 314 is configured to determine an approximation function that approximates the pricing parameter function 156 that yields substantially the same result as the original pricing parameter function. For example, the linearly transformed pricing parameter function can be expressed as:

${{{g(y)} \approx {q(y)}} = {\sum\limits_{j = 0}^{N_{E}}\; {\alpha_{j}{A_{j}(y)}}}},{y \in \left\lbrack {c,d} \right\rbrack}$

Where the collection {A_(j)(y)} is the desired family of orthogonal functions (321), the collection {α_(j)} is calculated by the approximation generator 314, and is given by:

α_(j) =<g(y)|A _(j)(y))

where (▪|▪) represents the scalar product that defines the orthogonality in the desired family of orthogonal functions. Typically, the {α_(j)} values are calculated via a numerical integration technique. For example, a Newton-Cotes methods, a trapezoidal method, a Simpson's method, a Boole's method, a Romberg integration method, a Gaussian quadrature methods, Chenshaw-Curtis methods, a Fejer method, a Gaus-Kronrod method, Fourier Transform methods, an Adaptive quadrature method, a Richardson extrapolation, a Monte Carlo and Quasi Monte Carlo method, a Markov chain Monte Carlo, a Metropolis Hastings algorithm, a Gibbs Sampling and Fast Fourier Transform methods.

In this case, the function q(y) is an accurate approximation of the linearly transformed pricing parameter function, g (y), when the upper bound of summation, NE, is sufficiently large. In mathematical notation:

ƒ(x)=g(y)≈q(y) for N _(E) sufficiently large

where y is the linear transform of x as per Step 364.

The approximation evaluator 316 can then determine the approximate value of the price parameter function at the desired time points and scenarios using the approximation function 319. The approximation evaluator 316 is used in place of the price parameter function evaluator 152 of the prior art forecasting system 100 when the price parameter scenarios are determined at Step 220 of the process carried out by the prior art forecasting system 100.

The approximation function 319 that approximates the pricing parameter function 156 is computationally simpler than the original pricing parameter function. Accordingly, the overall computational time required to carry out the step (Step 204) of determining the price parameter scenarios is reduced significantly. As discussed above, significant computational time needed in Step 204 in the prior art forecasting system 100 to achieve an accurate prediction, which can mean that the results of the prediction are no longer useful. As a consequence of determining the approximation function 319, the system 300 can significantly reduce the total computational time to generate the forecast, making the results useful.

The computational speed gained by the approximation evaluator 316 over the price parameter function evaluator 152 is particularly relevant for a number of families of orthogonal functions 319 in the approximating data function database 318, as for them either the gain in speed comes with no substantial loss of accuracy, or the system 300 can produce more accurate parameter values or metrics than 100 in the same time frame, or the system 300 can produce both an increase in speed and in accuracy over prior art.

An example of the process 348 carried out by the improved price parameter module 302 will now be described where the family of orthogonal functions can deliver an approximating function in the interval between a minimum value of −1 and a maximum value of 1 (i.e. yε[−1,1]), for example a Chebyshev polynomial.

Firstly, the transformer 312 transforms at Step 364 the pricing parameter function ƒ(x) which is in the interval [a,b] to an equivalent function in the interval [−1,1]. This is achieved starting with the equation of a straight line:

y=mx+c

Where m is the gradient of the straight line and c represents the intercept of the line with the y-axis. The mapping between the interval [a,b] to the interval [−1,1] means the line is known to pass through two points: the minimum value of x, (i.e. x_(min)=a), at y=−1 and the maximum value of x, (i.e. x_(max)=b), at y=1. Accordingly, it can be shown that:

$y = {{\frac{2}{b - a}x} - \frac{b + a}{b - a}}$

Expressing this linear transformation function as y=L(x), x can be expressed as the inverse of the function L such that x=L⁻¹ (y). In this example:

$x = {{L^{- 1}(y)} = \frac{{\left( {b - a} \right)y} + b + a}{2}}$

The transformer 312 can then determine the linearly transformed price parameter function such that:

ƒ(x)=ƒ(L ⁻¹(y))=g(y)

Once the transformer 312 has determined the linearly transformed price parameter function, the approximation generator 314 can generate at Step 380 a function that approximates the linearly transformed price parameter function based on a predetermined expansion series 321 from the expansion series database 320.

For example, one series that could be used to approximate the linearly transformed price parameter function, g(y), is:

${{c(y)} = {\sum\limits_{j = 0}^{N_{E}}\; {c_{j}{T_{j}(y)}}}},{y \in \left\lbrack {{- 1},1} \right\rbrack}$ T₀(y) = 1 T₁(y) = y T_(j + 1)(y) = 2yT_(j)(y) − T_(j − 1)(y), j ∈ ℕ  and  j > 1 ${c_{j} = {\frac{2}{\pi}{\int_{- 1}^{1}{\frac{1}{\sqrt{1 - y^{2}}}{g(y)}{T_{j}(y)}{y}}}}},{j > 0}$ ${c_{0} = {\frac{1}{\pi}{\int_{- 1}^{1}{\frac{1}{\sqrt{1 - y^{2}}}{g(y)}{y}}}}},{j = 0}$

The coefficients c_(j) can be computed by the approximation generator using a Fast Fourier Transform (FFT) function or any other numerical integration method known in the art. Typically, these methods have to call the original function (ƒ(x)) only in a relatively low number of times, in the “anchor” points in the ‘x’ domain.

In this example, the function c(y) is an accurate approximation of the linearly transformed pricing parameter function, g(y), for yε[−1,1], when the upper bound of summation, N_(E), is sufficiently large, and when g(y) is Lipschitz continuous in [−1,1]. It must be noted that the substantial majority of functions g(y) in practical applications have this Lipschitz property. In mathematical notation:

ƒ(x)=g(y)≈c(y) for N _(E) sufficiently large

Also, the approximation function (e.g., c(y)) can be calculated or approximated by an interpolation scheme, where the anchor points of it can be the zeros of each member of the chosen family of orthogonal functions, or a subset of the zeros.

Generally, in many forecasting systems, using this series expansion 321 a value of N_(E) (i.e. anchor points) between 1 and 200 is large enough to generate an approximation function, c(y), that outputs substantially the same result (typically the machine precision accuracy known in the art of 10⁻¹⁵) as the original price parameter function at the same time point and scenario. The accuracy is the difference between the output of the original price parameter function and the approximating function at the same set of input parameters. In some examples, the value of N_(E) may be in the range of 5 to 50 to give a better balance between computational time and accuracy of the results; in many cases this range provides an accuracy between 10⁻⁴ and 10⁻⁶, which is often considered as substantially zero error in the art. In some other cases, a relatively low N_(E) (e.g., 10) may give an accuracy of 10⁻², which can be good enough to produce meaningful metrics and subsequent forecasts in a time frame that makes them useful. By 10⁻², 10⁻⁴, 10⁻⁶, 10⁻¹⁵ or, generally, 10^(−Q) is meant that the lack of accuracy (e.g., the error) starts in the 2^(nd), 4^(th), 6^(th), 15^(th) or Q-th significant figure. The choice of N_(E) is a balance between the desire to minimise the error in the approximation function 319 and reduce the time taken to determine the desired number of price scenarios. In other words, choosing a larger value for N_(E) would minimise the error in the approximation function 319 but increase the time taken for the approximation generator 314 to generate an approximation function, and accordingly increase the overall time required to determine the scenarios. Conversely, choosing a smaller value for N_(E) would reduce the time taken to for the approximation generator 314 to generate an approximation function, and accordingly reduce the overall time require to determine the scenarios, however, the error in the output of the approximation function would increase. Clearly, when the number of scenarios is much greater than the N_(E) (for example, the number of scenarios, N, may be 20,000 whereas N_(E) may be 20), the process 348 carried out by the forecasting system 300 is much faster than the prior art process 200.

The approximation function 319 is then evaluated at Step 354 as many times as required to generate the scenarios as desired.

In addition to N_(E) and the number of scenarios (N_(S)), other factors that need to be accounted for to determine the improvement over prior art are the time it takes an approximation function, c(y), to run (t) and the time it takes the price parameter function to run (T). In general, the pricing step 204 of prior art systems will take approximately a total time N_(S)·T to run; the pricing step 354 in system 300 will take approximately N_(E)·T+N_(S)·t. System 300 will be beneficial over prior art when N_(S)·T is greater than N_(E)·T+N_(S)·t. That is the case in many forecasting systems. In other words, the value of N_(E) (i.e. the number of anchor points) that reduces the processing time required compared to the prior art system is:

$\frac{N_{S} \cdot T}{{N_{E} \cdot T} + {N_{S} \cdot t}} > 1$

Taking an example where the desired number of scenarios for a time point is ten thousand, where an N_(E) of 20 is enough to give a sufficiently good approximation, and where the time to compute the approximation function is lower than the price parameter function, in the prior art method, the price parameter function would be evaluated ten thousand times. However, using the approximation module 302, the approximation generator 314 would evaluate the linearly transformed price parameter function up to 20 times in order to determine an approximation function 319. The computationally faster approximation function 319 would be evaluated ten thousand times by the approximation evaluator 316 to determine the desired ten thousand price parameters for a given time step.

Choice of Orthogonal Function

The determination of which orthogonal function expansion series 321 (or approximations to them) from the expansion series database 320 is used to generate the approximation function 319 can be done in various ways. For example:

1. Practitioner experience

-   -   The user of the forecasting system 300 may predetermine the         desired orthogonal function series expansion 321 used by the         approximation generator 314 and makes this choice based on prior         knowledge and experience of using the forecasting system 300.         2. Trial and error     -   The user of the forecasting system 300 may run tests with a         plurality of the series expansions 321 and determine which         orthogonal function yields the most accurate result in the         desired computing time.         3. Run all and analyse     -   The approximation generator 314 of the forecasting system 300         may be configured to run all the series expansions 321 available         (or a subset of series expansions considered from historical         experience and data known to be the most likely to produce the         best results) to determine an estimate of the time required and         the associated error to complete the generation of a forecast.         The approximation generator 314 would then choose the optimum         orthogonal function for the forecast.

Improvements Achieved Over Prior Art

Discussed below are the results of the following tests comparing the prior art forecasting system 100 with the forecasting system 300 embodying the present invention:

1. Impact on the forecasting computing time for two standard pricers in a counterparty risk Monte Carlo simulations

-   -   a. Black-Scholes option pricer     -   b. Longstaff-Schwartz option pricer         2. Impact in computing time as a function of pricing function         computing time in a counterparty risk Monte Carlo simulation.         3. Reduction in noise in a counterparty risk Monte Carlo         simulation.

1. Computing Time

a. Black-Scholes Option Pricer

In an embodiment of the present invention, the p(y) series expansion was used to approximate the relatively simple Black-Scholes pricing parameter function in a forecasting Monte Carlo simulation, with 105 time points and 10,000 scenarios. In the computer where the embodiment of the invention was implemented, the time required to generate the 105×10000 price parameter scenarios using the price parameter module 106 of the prior art forecasting system 100 was 119.1 seconds, but the time was reduced to 25.6 seconds using the approximation module 302 of the forecasting system 300. This shows that the forecasting system 300 is 4.7 times faster than the prior art system to determine the pricing scenarios in Step 204.

It must be noted that the Black-Scholes pricing function used in the above test is already generally considered to be a super-fast function in the art. Accordingly, more common price parameter functions would result in a far greater speed-up multiple.

The overall time to carry out the process 200 and determine the counterparty credit risk metrics was 130.4 seconds using the prior art forecasting system 100 and 36.7 seconds using the forecasting system 300. This demonstrates a 3.6 times improvement in the overall time required to perform the same forecast.

FIGS. 11A and 11B show the EPE (Expected Positive Exposure, one of the counterparty credit risk metrics) profiles 400, 402 (one of the counterparty credit risk metrics) of this OTC derivative contract, both from the direct calculation of the price parameter function (profile 400) and via the approximation function (profile 402). The figures show the EPE profiles for different collateralisation conditions. It can be seen that the results for both collateralised and uncollateralised cases are substantially the same (within the always existent numerical noise of the Monte Carlo simulation).

b. Longstaff-Schwartz Option Pricer

This is a typical pricer function used for American options. These options are fairly common, but pricing is too slow using the prior art forecasting system 100 for good counterparty risk calculations.

For example, in the computer in which this tests were carried out, the time to compute a single Black-Scholes option price was around 0.6 milliseconds, whilst that time for a Longstaff-Schwartz option in this test was around 130 milliseconds. Again, the test was performed with 105 time points and 10,000 scenarios. The time to compute an approximated single price using the present embodiment in these tests was around 0.15 milliseconds.

The time required to generate the price parameter scenarios using the prior art forecasting system 100 was 25,112 seconds but was reduced to 57.5 seconds using the forecasting system 300 comprising the approximation module 302. This shows that the forecasting system 300 of the present embodiment is 437 times faster than the prior art forecasting system 100 to determine the pricing scenarios in Step 204. The overall time to carry out the process 200 and determine the counterparty credit risk metrics was 25,134 seconds in the prior art forecasting system 100 and 69.3 seconds using the forecasting system 300 of the present embodiment. This demonstrates a 363 times improvement in the overall time required to perform the same calculation.

FIGS. 12A and 12B shows the PFE profiles (Potential Future Exposure, one of the counterparty credit risk metrics) of this OTC derivative contract, both from the direct calculation of the price parameter function (profile 404) and via the approximation function 319 (profile 406). Two plots are provided for the collateralised and uncollateralised cases. It can be seen that the results for both the collateralised and uncollateralised cases are substantially the same (within the always existent numerical noise of the Monte Carlo simulation).

2. Impact as a Function of Pricing Time

This performance test involves using a simple generic pricing function as possible comprising a single line of code with a simple subtraction. A time delay was introduced into the original function. The time delay was varied and the impact on the computing time required is shown in the table below, wherein the tests were carried out under the same conditions using the same computer hardware.

Single price Price parameter scenarios determination Overall risk metric determination parameter time time function Prior art Approximation Prior art Approximation time system system system system (millisecond) (seconds) (seconds) Improvement (seconds) (seconds) Improvement 0.03 6.5 28.6 0.2 20.4 42.2 0.5 1.03 198 28.5 6.9 211 41.8 5.0 5.03 952 28.5 33 966 41.8 23 10.03 1903 29.2 65 1917 42.7 45 50.03 9511 30.0 317 9525 43.5 219 100.03 19021 29.7 641 19035 43.1 442

This data is graphically shown in FIG. 13 where the prior art direct pricing graph 408 shows a relationship of linearly increasing computing time with linearly increasing price parameter function computing time. Conversely, the pricing graph 410 generated by the present embodiment shows an almost constant simulation computing time with increasing individual price parameter function computing time. Financial institutions can easily have OTC derivatives that have a pricing time in the excess of 1 second; in some cases of several minutes, or even a few hours in some extreme cases. Extrapolating the results data in the table above, with a function pricing time of 1000 milliseconds, for example, the forecasting system 300 could improve the computing time required by 5000 times.

Also, in the prior art, it is useless to try to calculate forecasts from counterparty credit risk metrics for an OTC derivative that has a pricing function time of 1 minute. However, with embodiments of the present invention, the pricing calculation is achievable in a very reasonable time frame.

3. Reduction in Noise

This test uses the Longstaff-Schwartz pricing function where 100 scenarios are run using the prior art forecasting system 100 and 400,000 scenarios are run using the forecasting system 300. Both tests used 105 time steps.

The time taken for the process 204 with direct evaluation of the pricing parameter function (namely the prior art method) was 13,955 seconds. The time taken using the approximation approach for the forecasting system 300 was comparable at 12,248 seconds. However the noise in the PFE results 412 from this simulation with direct evaluation of the pricing parameter function of the prior art was noticeably higher than the noise 414 generated by the forecasting system 300 of the present embodiment, as it can be seen in FIGS. 14A and 14B.

This test shows that, for the approximately the same computational time, a more accurate and less noisy risk metric can be determined using the forecasting system 300 over the prior art forecasting system 100.

Application in Multiple Dimensions

The embodiments described above have been discussed with reference to a one dimensional price parameter function, i.e. with a single variable, x. In other embodiments, the price parameter function 156 may be multidimensional i.e. multivariable, having n-dimensions, e.g. ƒ(x₁,x₂,x₃, . . . , x_(n)). The forecasting system 300 can handle multidimensional price parameter functions in a number of ways. Four methods that may be carried out by the forecasting system 300 to perform multi-dimensional forecasting are discussed below by way of example only. It is to be understood that the forecasting system 300 may be configured to carry out multi-dimensional forecasting in one or more of several ways:

1. Locking the Price Parameter Function to One Degree of Freedom

In most embodiments, one of the variables out of x₁ to x_(n) causes most of the change in the output of a price parameter function 156. In this case, the value of the other variables can be kept constant as changing the value of the variable generally has a negligible effect on the output of the price parameter function 156. Using the notation x_(i)={tilde over (x)}_(i), to represent a variable being held constant, taking the example where x₁ is the dominant variable in a price parameter function 156, it can be stated that:

ƒ(x ₂ ,x ₃ , . . . x _(n))≈ƒ(x ₁ ,{tilde over (x)} ₂ ,{tilde over (x)} ₂ . . . x _(n))

As can be seen, the number of degrees of freedom locked is one less than the number of dimensions. Using this approximation, the approximation generator 314 could determine the approximating function 319 in the one-dimensional approach of Step 354 as described above with reference to FIGS. 6A and 6B.

For example, taking a price parameter function that outputs the temperature in a weather forecasting system, where the price parameter function comprises two variables, humidity (x₁) and height (x₂), when forecast the temperature over an area of land that is fairly flat, the height, x₂, is taken to be a constant value, {tilde over (x)}₂. The function (x₁,{tilde over (x)}₂) is then approximated by the forecasting system 300.

2. Coordinate Transformation

In other embodiments, the coordinate system is changed by the forecasting system 300 such that the method described above in locking the price parameter function to one degree of freedom can be applied. For example, taking a two variable case where x₁+x₂ causes the predominant change in the output of the price parameter function, the following change of coordinates would be applied:

z ₁ =x ₁ +x ₂

z ₂ =x ₁ −x ₂

Such that, in the new frame of reference:

ƒ(z ₁ ,z ₂)≈ƒ(z ₁ ,{tilde over (z)} ₂)

This price parameter function can then be approximated by the approximation generator 314 as described above.

Whilst this example presents a two-dimensional price parameter function, it is to be understood that this can be readily applied to a general n-dimensional price parameter function without loss of generality.

3. Approximating in More than One Dimension

In some embodiments, a plurality of the variables may have a strong influence on the output of the price parameter function 156 in a combination that makes it very computationally difficult to apply the coordinate transformation and/or locking the price parameter function to one degree of freedom in the methods described above.

Taking an example two-dimensional price parameter function, ƒ(x₁, x₂), where the variable x₁ causes the majority of the changes in the output of the price parameter function, but the variable x₂ also causes non-negligible changes in the output of ƒ(x₁, x₂). By taking two representative values of x₂, e.g. the minimum and maximum values ({tilde over (x)}₂ ^(min), {tilde over (x)}₂ ^(max)), two one-dimensional approximation functions can be solved locking the price parameter function to one degree of freedom method explained above for ƒ_(min)(x₁, {tilde over (x)}₂ ^(min)) and ƒ_(max)(x₁, {tilde over (x)}₂ ^(min)). Then a linear interpolation in the variable x₂ is used to find out a good approximation for ƒ(x₁,x₂).

A fully general way of solving this problem would be by transforming the n-dimensional approximation problem to n one-dimensional approximation problems.

A two-dimensional embodiment is explained as an example below, however, it is to be understood that the same method can readily be extended to n-dimensions.

Starting from the one-dimensional example of a series expansion described above:

${{c(y)} = {\sum\limits_{j = 0}^{N_{E}}\; {c_{j}{T_{j}(y)}}}},{y \in \left\lbrack {{- 1},1} \right\rbrack}$ T₀(y) = 1 T₁(y) = y T_(j + 1)(y) = 2yT_(j)(y) − T_(j − 1)(y), j ∈ ℕ  and  j > 1 ${c_{j} = {\frac{2}{\pi}{\int_{- 1}^{1}{\frac{1}{\sqrt{1 - y^{2}}}{g(y)}{T_{j}(y)}{y}}}}},{j > 0}$ ${c_{0} = {\frac{1}{\pi}{\int_{- 1}^{1}{\frac{1}{\sqrt{1 - y^{2}}}{g(y)}{y}}}}},{j = 0}$

Firstly, a linear transformation is performed on each variable, x₁ and x₂, by the transformer 312 of the approximation module 302, such that:

x ₁ ε[a ₁ ,b ₁ ]→y ₁ ε[c ₁ ,d ₁]

x ₂ ε[a ₂ ,b ₂ ]→y ₂ ε[c ₂ ,d ₂]

Where, in this example, c₁=c₂=−1, d₁=d₂=1.

The two one-dimensional approximations are then built each with an approximation order of N_(E4) and N_(E,2). FIG. 15 shows a graph 500 where N_(E,1)=3 and N_(E,2)=4, hence there are four primary anchor points 502 for variable y₁, and five primary anchor points 502 for variable y₂ for the calculation of the orthogonal approximation interpolant function. The anchor points are illustrated by the solid circles in FIG. 15.

The output of the price parameter function at the point X (labelled as 504 in FIG. 15), given by the values ({tilde over (x)}₁,{tilde over (x)}₂)→({tilde over (y)}_(x,1),{tilde over (y)}_(x,2)) can then be determined in the following steps:

-   -   1. Generate the information for the approximated evaluation     -   This is done by the approximation generator 314 which evaluates         ƒ(x₁,x₂)=g(y₁,y₂) at each primary anchor point 502. This         evaluation is done solving the full price parameter function,         and in the example of FIG. 15, it is evaluated 20 times.     -   2. Compute ƒ({tilde over (x)}₁,{tilde over (x)}₂)=g({tilde over         (y)}_(x,1),{tilde over (y)}_(x,2))     -   a) If ({tilde over (y)}₁,{tilde over (y)}₂) refers to the         primary anchor points 502, using them, an approximated value for         g({tilde over (y)}_(x,1),{tilde over (y)}₂) can be determined by         the approximation evaluator 316 to determine a set of five         secondary anchor points 506, illustrated by the hollow circles         in FIG. 15. One degree of freedom of the price parameter         function has been locked in order to determine the secondary         anchor points 506. In other words, a one-dimensional         approximation method has been determined for each secondary         anchor point 506, along y₁, in which the value of y₂ is held         constant at the value of each of the primary anchor points of         {tilde over (y)}₂, respectively.     -   b) Using the secondary anchor points 506, an approximation         function can be used by the approximation evaluator 316 to         determine the output of g({tilde over (y)}_(x,1),{tilde over         (y)}_(x,2)) along y₂, in which the value of y₁ is held constant         at the value of {tilde over (y)}_(x,1).     -   Notably, the determination of the secondary anchor points 506 in         the step of computing ƒ({tilde over (x)}₁,{tilde over         (x)}₂)=g({tilde over (y)}_(x,1),{tilde over (y)}_(x,2)) is fast         as it uses an approximation function 319, the advantages of         which have been described above.         4. Combination of Approximating in More than One Dimension and         Locking Degrees of Freedom

In other embodiments, the price parameter function 156 comprises a plurality of variables in which some variables have hardly any impact on the output of the price parameter function, but other variables do. In this case, a combination of both locking one or more degrees of freedom in the price parameter function and approximating in more than one dimension as described above can be used. Firstly, the non-negligible variables are not locked:

ƒ(x ₂ ,x ₃ , . . . x _(n))≈ƒ(x ₁ ,x ₂ , . . . x _(s) ,{tilde over (x)} _(s+1) , . . . ,{tilde over (x)} _(n))

Taking the number of non-negligible variables to be s, the number of degrees of freedom locked is s less than the number of dimensions. Then, the multi-dimensional technique can be applied to those variables that are not locked. Improvements Achieved Over Prior Art of Approximating in More than One Dimension

Taking an example where there are 105 time steps and 10,000 scenarios and taking the computation time to be 130 milliseconds (as in the case of the Longstaff-Schwartz example discussed above) for each evaluation of a two-dimensional price parameter function, the time taken to evaluate the price parameter function 156 for all time points and scenarios is:

105×10,000×130 milliseconds 37.9 hours

However, using the multi-dimensional technique in the two-dimensional price parameter function of an embodiment of the present invention, and assuming that N_(E,1)=N_(E,2)=10 (leading to 121 anchor points 502), and assuming the estimated computational time of each one-dimensional approximation function of 0.15 milliseconds (like in the Longstaff-Schwartz example), the estimated time to compute the same number of time points and scenarios is given by the time taken to evaluate price parameter function for anchor points added to the time taken to evaluate the secondary anchor points 506 for all time points and scenarios is:

11×11×130 milliseconds+105×10,000×(11+1)×0.15 milliseconds=0.5 hours

This yields an improvement of around 76 times less computing time.

Many modifications may be made to the above examples without departing from the spirit and scope of the present invention as defined in the accompanying claims.

For example, the price parameter function could be a function that can be used to evaluate the volume of crude oil extractable from an oil reservoir at a particular time in the oil and gas industry, or it can be used to determine the temperature in a location at a particular time in the field of weather modelling.

The functionality and processes described above may be implemented via a system architecture such as that shown in FIGS. 1A, 1B, 2A, 2B, 4, 5A, 5B, 6A, 6B and 6C. At its most general level, the forecasting system 300 according to the present invention comprises the risk factor module 104, the risk metric module 108, the approximation module 302 and the controller 102 which are each implemented by a processor or can be software instances implemented on any parallel processing processor platform. Example processing architecture that may be suitable for the present invention include the Intel Xeon, Core i7 and Core i5 x86-64 microprocessors utilising the Haswell, Sandy Bridge or Ivy Bridge microarchitectures, the AMD Opteron x86-64 microprocessors utilising the Bulldozer, Piledriver or Jaguar microarchitectures, or 64-bit ARM System-on-Chip architecture.

The forecasting system 300 can equally be used in different forecasting applications, other that OTC derivatives in financial modelling. For example, the forecasting system 300 could be adapted for use in molecular biology, oil and gas, structural engineering, and weather forecasting by including appropriate parameter functions into the pricing function database 154. Additionally, appropriate stochastic input variables would be determined by the risk factor module 104 and stored in the stochastic model database 112.

More specifically, returning to the example of oil and gas described in the background section, the set of one or more input parameters may include: viscosity of the crude oil, remaining supplies in the oil reservoir, oil pumping capacity, well temperature, availability and/or effectiveness of drilling fluid which each vary stochastically over time. In this example, the parameter database comprises a parameter function to which outputs the extractable volume of crude oil in an oil well, the parameter function using the input parameters. The processes 348 and 354 carried out by the forecasting system 300 are substantially the same compared to using a price parameter function. However, at Step 380, an approximation function is generated for the parameter function for outputting the extractable volume of crude oil in an oil well. The approximation function is less computationally expensive to evaluate relative to the original parameter function. The approximation function is evaluated for each time point in each scenario. The scenarios determined at Step 354 using the approximation function can then be used to determine metrics which may include required oil refinement and purification capacity and staffing levels required in downstream processes. Preferably, the approximation function used is based on fitting a set of orthogonal functions or an approximation to a set of orthogonal functions to the parameters.

Fitting, (i.e. curve fitting) used to generate the approximation function is the process of constructing a curve, or mathematical function, that has the best fit to a series of data points, possibly subject to constraints.

Accuracy of the approximation function or of a metric is the degree to which the result of a measurement, calculation, or specification conforms to the correct value or a standard.

Similarly, for weather forecasting, the set of one or more input variables may include air pressure, air density, solar energy and humidity each from hundreds or more of different locations. The parameter database comprises a parameter function to output the temperature at a geographic location as a function of one or more of the input variables. The processes 348 and 354 carried out by the forecasting system 300 are substantially the same compared to using a price parameter function. However, at Step 380, an approximation function is generated for the parameter function for outputting the temperature at a geographic location. The approximation function is less computationally expensive to evaluate relative to the original parameter function. The approximation function is evaluated for each time point in each scenario. The scenarios determined at Step 354 using the approximation function can then be used to determine metrics which may include sea level, pollen count, UV index and wind speed/direction. Preferably, the approximation function used is based on fitting a set of orthogonal functions or an approximation to a set of orthogonal functions to the parameters.

User devices may further be in communication with the forecasting system 300, e.g. via a communications network such as the Internet, and interface with the input/output module via a web-enabled browser or application running on a portable multifunctional device (e.g., smart phone, personal digital assistant (PDA), iPad™, notebook, subnotebook, tablet PC, electronic reader, etc.) or a non-portable device (e.g., desktop computer, such as PC or Mac, internet capable electronic displays/TVs, etc.).

The forecasting system 300 may further comprise one or more processor modules. In FIGS. 1A and 5A, five modules are shown. However it is to be appreciated that other arrangements of processing modules and databases may be envisaged which distribute the functionality of the system which would also be within the scope of the present invention.

Computer programs based on the written description and flow charts of this invention are within the skill of an experienced developer. The various software components, functions, programs and/or program modules described in this application may be created using any of the techniques known to one skilled in the art or can be designed in connection with existing software, and implemented as software code to be executed by a processor using any suitable computer language, such as, for example, Java®, C++, Perl, HTML or XML. The software code may be stored as a series of instructions or commands on a computer readable medium, such as a random access memory (RAM), a read only memory (ROM), a magnetic medium such as a hard-drive or a floppy disk, or an optical medium such as a CD-ROM. Any such computer readable medium may reside on or within a single computational apparatus, and may be present on or within different computational apparatuses within a system or network. 

1. A computer-implemented method of determining an approximated value of a parameter in a first domain, the parameter being dependent on one or more variables which vary in a second domain, and the parameter being determined by a function which relates sets of values of the one or more variables in the second domain to corresponding values in the first domain, the method being implemented on a computer system including a processor, and the method comprising: determining a plurality of anchor points in the second domain, wherein each anchor point comprises a set of values of the one or more variables in the second domain; evaluating, at each anchor point, the function to generate corresponding values of the parameter in the first domain; generating an approximation function to the function by fitting a series of orthogonal functions or an approximation to a series of orthogonal functions to the corresponding values of the parameter in the first domain; and using the approximation function to generate the approximated value of the parameter in the first domain.
 2. The computer-implemented method of claim 1, further comprising: transforming the function into an intermediary function before the generating step, wherein the intermediary function is a linear transformation of the function, and wherein the generating step comprises approximating the intermediary function.
 3. The computer-implemented method of claim 2, wherein the transforming step comprises linearly transforming the function into a predetermined range.
 4. The computer-implemented method of claim 1, wherein the series of orthogonal functions is one of the group comprising: sines, cosines, sines and cosines, Bessel functions, Gegenbauer polynomials, Hermite polynomials, Laguerre polynomials, Chebyshev polynomials, Jacobi polynomials, Spherical harmonics, Walsh functions, Legendre polynomials, Zernike polynomials, Wilson polynomials, Meixner-Pollaczek polynomials, continuous Hahn polynomials, continuous dual Hahn polynomials, a classical polynomials described by the Askey scheme, Askey-Wilson polynomials, Racah polynomials, dual Hahn polynomials, Meixner polynomials, piecewise constant interpolants, linear interpolants, polynomial interpolants, gaussian process based interpolants, spline interpolants, barycentric interpolants, Krawtchouk polynomials, Charlier polynomials, sieved ultraspherical polynomials, sieved Jacobi polynomials, sieved Pollaczek polynomials, Rational interpolants, Trigonometric interpolants, Hermite interpolants, Cubic interpolants, and Rogers-Szegö polynomials.
 5. The computer-implemented method of claim 1, wherein the series of orthogonal functions is a series of orthogonal polynomials.
 6. The computer-implemented method of claim 1, wherein the approximation to a series of orthogonal functions is an approximation to a series of orthogonal polynomials.
 7. The computer-implemented method of claim 1, further comprising: selecting a series of orthogonal functions or an approximation to a series of orthogonal functions before the generating step, and wherein the generating step comprises using the selected series of orthogonal functions or the approximation to the series of orthogonal functions.
 8. The computer-implemented method of claim 1, wherein the number of anchor points, N_(E), is so that: $\frac{N_{S} \cdot T}{{N_{E} \cdot T} + {N_{S} \cdot t}} > 1$ Wherein N_(S) is the number of scenarios, t is the time taken to run the approximation function, and T is the time taken to run the function.
 9. The computer-implemented method of claim 1, further comprising: holding the values of all but one variable in the set to be constant in the evaluating step.
 10. The computer-implemented method of claim 1, further comprising using the approximated values of the parameter to generate standard metrics.
 11. The computer-implemented method of claim 1, wherein the values of the variables vary stochastically.
 12. The computer-implemented method of claim 1, wherein the using step comprises using the approximation function a plurality of times to generate at least one scenario or one time step, each scenario or each time step comprising a plurality of approximated values.
 13. The computer-implemented method of claim 12, wherein the number of scenarios or time steps is significantly greater than the number of anchor points.
 14. The computer-implemented method of claim 1, wherein an output of the function is a parameter of a financial product.
 15. The computer-implemented method of claim 14, wherein the financial product is a financial derivative including one or more of a group comprising: an option pricing function, a swap pricing function and a combination thereof.
 16. The computer-implemented method of claim 1, wherein the function is one of the group comprising: a Black-Scholes model, a Longstaff-Schwartz model, a binomial options pricing model, a Black model, a Garman-Kohlhagen model, a Vanna-Volga method, a Chen model, a Merton's model, a Vasicek model, a Rendleman-Bartter model, a Cox-Ingersoll-Ross model, a Ho-Lee model, a Hull-White model, a Black-Derman-Toy model, a Black-Karasinski model, a Heston model, a Monte Carlo based pricing model, a binomial pricing model, a trinomial pricing model, a tree based pricing model, a finite-difference based pricing model, a Heath-Jarrow-Morton model, a variance gamma model, a Fuzzy pay-off method, a Single-index model, a Chepakovich valuation model, a Markov switching multifractal, a Datar Mathews method, and a Kalotay-Williams-Fabozzi model.
 17. The computer-implemented method of claim 1, wherein the anchor points are the zeros of each of orthogonal function, or a subset of the zeros.
 18. The computer-implemented method of claim 17, wherein the approximating function is generated using an interpolation scheme.
 19. The computer-implemented method of claim 18, wherein the interpolation scheme is one from the group comprising: piecewise constant interpolants, linear interpolants, polynomial interpolants, gaussian process based interpolants, spline interpolants, barycentric interpolants, Rational interpolants, Trigonometric interpolants, Hermite interpolants and Cubic interpolants.
 20. The computer-implemented method of claim 1, wherein the anchor points are integration points of a numerical integration scheme.
 21. The computer-implemented method of claim 20, wherein the numerical integration scheme is one from the group comprising: Newton-Cotes methods, a trapezoidal method, a Simpson's method, a Boole's method, a Romberg integration method, Gaussian quadrature methods, Chenshaw-Curtis methods, a Fejer method, a Gaus-Kronrod method, Fourier Transform methods, an Adaptive quadrature method, a Richardson extrapolation, a Monte Carlo and Quasi Monte Carlo method, a Markov chain Monte Carlo, a Metropolis Hastings algorithm, a Gibbs Sampling, and Fast Fourier Transform methods.
 22. The computer-implemented method of claim 1, wherein the speed of the calculation is increased with no loss of accuracy in the standard metrics, and/or the accuracy of the standard metrics increase with no decrease in the speed of the calculation when using the approximation function.
 23. The computer-implemented method of claim 1, wherein there is a difference between the approximated value of the parameter and the parameter in the 2^(nd) significant figure when compared at the same point in the first domain.
 24. The computer-implemented method of claim 1, wherein there is a difference between the approximated value of the parameter and the parameter in between the 4th and the 6^(th) significant figure when compared at the same point in the first domain.
 25. The computer-implemented method of claim 1, wherein there is a difference between the approximated value of the parameter and the parameter in the 15^(th) significant figure when compared at the same point in the first domain.
 26. A financial derivative comprising a parameter, wherein a value of the parameter is determined using the computer-implemented method of claim
 1. 27. A computer system comprising a processor configured to determine an approximated value of a parameter in a first domain, the parameter being dependent on one or more variables which vary in a second domain, and the parameter being determined by a function which relates sets of values of the one or more variables in the second domain to corresponding values in the first domain, the computer system comprising: a determination module arranged to determine a plurality of anchor points in the second domain, wherein each anchor point comprises a set of values of the one or more variables in the second domain; an evaluation module arranged to evaluate, at each anchor point, the function to generate corresponding values of the parameter in the first domain; a generation module arranged to generate an approximation function to the function by fitting a series of orthogonal functions or an approximation to a series of orthogonal functions to the corresponding values of the parameter in the first domain; and an approximation module arranged to use the approximation function to generate the approximated value of the parameter in the first domain.
 28. A computer-implemented method of determining an approximated value of a parameter in a first domain, the parameter being dependent on one or more variables which vary in a second domain, and the parameter being determined by a function which relates sets of values of the one or more variables in the second domain to corresponding values in the first domain, the method using a series of orthogonal functions or an approximation to a series of orthogonal functions to approximate the function. 